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The calculation of the self force in the modeling of the gravitational- wave emission from extreme- 
mass-ratio binaries is a challenging task. Here we address the question of the possible emergence 
of a persistent spurious solution in time-domain schemes, referred to as a Jost junk solution in 
the literature, that may contaminate self force calculations. Previous studies suggested that Jost 
solutions are due to the use of zero initial data, which is inconsistent with the singular sources 
associated with the small object, described as a point mass. However, in this work we show that 
the specific origin is an inconsistency in the translation of the singular sources into jump conditions. 
More importantly, we identify the correct implementation of the sources at late times as the sufficient 
condition guaranteeing the absence of Jost junk solutions. 

PACS numbers: 04.25.dg, 04.30.Db, 95.30.Sf 



Extreme-mass-ratio inspirals are one of the most im- 
portant sources of gravitational radiation for the future 
space-based observatory Laser Interferometer Space An- 
tenna (LISA). They consist of a stellar compact object 
inspiralling into a massive black hole (MBH; with mass 
in the range M = 10^ — IQ^ Mq). They are long lasting 
sources that in the last year before plunge can spend 
of the order of 10^ cycles inside the LISA frequency 
band [1]. To extract these signals from the future LISA 
data stream we require very precise theoretical waveform 
templates, as the signal-to-noise ratio accumulates slowly 
with time [2]. To achieve such precision we need a very 
accurate description of the slow inspiral, which can be 
seen as due to the action of a local force, the self force. 
The computation of the self force is a challenge since we 
need to regularize the gravitational perturbations created 
by the stellar compact object, which is modeled as a point 
mass. Nevertheless, during the last decade there has been 
significant progress in the self-force program (see for 
reviews) and different methods to compute it have been 
developed. Of particular relevance is the computation of 
the gravitational self force in the case of a nonrotating 
MBH, both for circular and eccentric orbits OEj, using 
time-domain methods. 

The aim of this letter is to clarify an issue related to 
time-domain self-force calculations: The appearance of 
the so-called Jost junk solutions discussed recently in [7] , 
which can contaminate time-domain calculations of the 
self force, limiting their accuracy. Here, following [7j, we 
show that such solutions are not a generic feature of time- 
domain schemes, but rather a result of the use of trivial 
zero initial data combined with a particular implemen- 
tation of the singular sources associated with the point 
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mass. We then identify the specific ingredients responsi- 
ble for the eventual appearance of the spurious solution 
and point out a straightforward solution, eliminating in 
this way concerns about the generation of this kind of 
persistent junk contamination when calculating the self 
force. 

The self force is determined by the gravitational per- 
turbations created by the point mass as it orbits the 
MBH, what we call the retarded field, whose computa- 
tion is the main numerical task in time-domain self-force 
schemes. The retarded field is singular at the location of 
the point mass and needs to be regularized by using ana- 
lytic expressions for the singular field (see, e.g. ^Sj). The 
retarded field is described by a system of coupled wave- 
type equations that derive from the perturbed Einstein's 
equations, which means that their structure depends on 
the gauge adopted. Not all the gauges are suitable for 
self-force computations. 

The principal part of the wave- type equations that gov- 
ern the retarded field is common to most gauges. For 
nonrotating MBHs, each spherical harmonic component 
obeys decoupled 1-1-1 wave equations whose structure is 
captured by the following model equation: 

[-d^^ + dl. -V{r)\^{t,r) = S{t,r) 

= /(r) [ G{t, r) 6{r - r^{t)) + F(i, r) <5'(r - r^it)) ](1) 

where /(r) = 1 - 2M/r, M is the MBH mass, r* is the 
tortoise coordinate r* = r + 2Mln[r/(2M) — 1], and rp(i) 
is the particle's radial motion. This model equation con- 
sists of a 1 -I- 1 wave operator —d^ -f 9^, , a potential term, 
F(r)4', and a singular source term S{t,r). Equation ([T]) 
has the form of the master equations that rule perturba- 
tions of different fields in Schwarzschild spacetime. They 
are tied to a particular set of gauges, but in other gauges 
the principal structure of the equations that govern the 
dynamics of the perturbations is similar. Therefore, we 
will use the model Eq. ([T]) in order to analyze the appear- 
ance of the Jost junk solution, as done in [71. Without 
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loss of generality, we deal with the case of circular orbits 
(j-p = const). 

Time-domain schemes address the resolution of Eq. ([T]) 
as an initial-boundary value problem, by prescribing ini- 
tial conditions on a time slice together with appropriate 
boundary conditions along the evolution. Appropriate 
initial data (ID) at a finite time are not known for Eq. ([T]). 
A common practice consists of setting trivial ID, namely 



(2) 



This ID is inconsistent with the singular structure of the 
source and corresponds to a solution that is not contin- 
uous in time at t — to. As a consequence, when solving 
(numerically) Eq. ^ an initial burst of junk radiation is 
produced and one must wait until it has been radiated 
away in order to compute the self force. This strategy 
relies on the assumption that junk radiation is actually 
radiated away. In this sense, Ref. [7] addresses the ques- 
tion of whether the use of trivial ID can give rise to spu- 
rious solutions persistent in time. To answer this ques- 
tion a double complementary numerical and analytical 
approach is adopted in [7] that we reanalyze below. 

(a) Numerical approach. The impact of inconsistent 
ID is assessed by constructing, first, a solution to Eq. 
([T| with trivial ID, referred to as Vfimpuisivo- Second, the 
sources m Eq. are modified to make them compatible 
with trivial ID. This is achieved by a smooth switch on 
in time of the sources 

F^it, r) EE T)F[t, r) , G^(i, r) = a{t, T)G{t, r) , (3) 

where a(i,T) smoothly interpolates between and 1, at 
initial and late times, i.e. 0(^0,^) = and a(t, r) — 1 
(for t ^ r). The solution obtained with trivial ID and 
smooth sources ([s]) is referred to as 'I'sinooth- At late 
times, f ^ T, a numerical function ^'^Qgj. is defined as 



^Jost = ^Impulsive — ^Smooth • 



(4) 



The function ^^^^^ has the following properties [HI: (i) 
It is time independent: c^t^j^^gj = 0, Vi. (ii) It has a jump 
at the particle: [*jost]p = ~fp^F{to,rp), where fp = 

f{rp). (iii) The spatial derivative, 9r»^'j'ost, is continuous 
at r = Tp. 

(b) Analytical approach. Motivated by the numerical 
approach, one constructs the analytical function ^^^^t ^ 



*Jost ^ ^Jost ^+ 1 



(5) 



where 6+ = Q{r* - r*p), e_ = e(r; - r*), being 6 the 

Heaviside step function, and ^^^^ and ^'^o^j solve the 
homogeneous, stationary version of Eq. ([ij 



[d% - y(r)] * 



A,± 
Jost 



(6) 



such that [^'j'ostlp = ^-fp ^F{to,rp). This means (see 
below) that ^^^gt satisfies an inhomogeneous version 
of Eq. ([g]) with a stationary singular source given by 



Differences between vE'j'ost and 



-fpF{to,rp)5'{r - rp) 
^'j'ggj are shown to vanish within numerical precision 
in [7]. Thus, we refer to a single Jost function, 5'jost- 

From this analysis we conclude that at late times, 
when the time switch-on function a{t,T) equals 1 and 
the sources for 5'impuisivc and 5'smooth coincide, their dif- 
ference should be a solution to the homogeneous version 
of Eq. ([iJ. However, this is in conflict with the fact that, 
as discussed above, ^'jost solves a (stationary) version 
of Eq. (Tj) with a singular distributional source. Cer- 
tainly, this contradiction arises from the use of inconsis- 
tent sources and ID. But more importantly, it also sug- 
gests strongly that, at late times, one is actually solving 
two different systems for ^'impulsive and ^'smooth, namely 
with different sources. To assess this point we now dis- 
cuss the implementation of Eq. ([l]) in jTJ . 

In the particle-without-particle (PwP) approach to 
Eq. ( [I]), introduced in [9l [10] and further discussed 
in [Till 11 [T^ . the point mass is placed at the bound- 
ary between two integration domains. Then, one solves 
a homogeneous problem in each domain and the singu- 
lar sources in Eq. are translated into jump condi- 
tions at the particle. To illustrate this it is convenient 
to rewrite Eq. ([T]) as a first-order hyperbolic system (cf. 
Refs. [121 US] for a second-order formulation) . To that 
end, we introduce the fields [TS] 

(f> = dt^ , ip = dr*'^ , (7) 

and Eq. becomes a system for the vector (^f, (/), ip) 

dt-9 = (b, dt(j}^dr'ip-Vir)^-Sit,r), dtv^dr*dp.{%) 

In the PwP approach we perform the splitting 



+ ' 



(9) 
(10) 



^ = (^-e_ + ^+e+ + [vi']^5(r*-r;), (ii) 

where the Dirac term of ip follows from consistency be- 
tween the definition of tp in ( [7| a nd the expression for 
^ in ([9]). Inserting expressions |9|-(11) in system ([s]) we 
obtain homogeneous systems of equations for the fields 



dr'P^ - V(r)-^^ , dtp^^ = 9,.0±(,12) 



and jump conditions to communicate them across the 
particle 



[*]pW 



fn'F{t,rp), 



ip 

p-i 



[4>]p{t) = f-'id,F)it,rp), 
[^]^{t) = G{t,rp)~fp\d,,F){t,rp). 



(13) 
(14) 
(15) 



This whole system, i.e. homogeneous Eqs. (12 1 and 
jumps (13l-(15), is equivalent to the original Eq. dl]). 
Two remarks are in order, regarding the system ( |12p - 
(15). First, not all initial conditions are consistent with 



the source S{t, r), since trivial ID violates the jump con- 
ditions (13)-(15|. Second, jumps [^]p and [0]^ contain 
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redundant information since, consistently with the defini- Then, the equations for (0, (p) become 



tion of in (l7|), is the time derivative of [ 5* ] . How- 



ever, it is crucial to realize that the jump condition (14) 
for does not account for the initial value of ['5]^. 
Therefore, if the evolution scheme fails to implement the 
condition 



(16) 



it is no longer equivalent to the model Eq. ([T]) . In [71 [T^] , 

this equation is written in first-order form in terms of the 



variables (^E*, 0, 1^9), as in (13l-(15), but the equations for 
<f> and are now 



5f0 dr'ip~V{r)-^ - J^5{r* 
dtip = dr'(t>+ J^5{r* ~ r*) , 



(17) 



where J<p = [v]p ^nd — [<t']p- The jump conditions 
can be recovered by multiplying the equations by a test 
function, integrating (by parts) between r* — e and r* +e, 
and taking the limit e 0. However, system (17 1 does 
not explicitly enforce the jump condition on [ ^' ]^ More 
precisely, the jump condition on [ ^ ]p is expected to be 



implemented by enforcing 



which means that [ ^ 1 



is enforced up to an initial condition: Only the value of 
(9t [ ^' ]p is enforced. Since 5* is coupled to the rest of 
the system through the potential term, V{r)'i>, this has 
consequences on the whole system. In conclusion, the 
evolution system (17) is not equivalent to Eq. ([!]). 



Another observation is that trivial ID is confusing: 
From the source perspective we must impose the con- 



dition ( 16 ), but from the point of view of the ID it seems 
reasonable to choose [ 5" ]p (to) = 0. This confusion is just 
a consequence of the inconsistency between the singular 
source term and trivial ID that we pointed out above. 

The main observation of this work is that Jost junk 
solutions appear as a consequence of implementing a 
finite jump condition, [^]p, by enforcing an infinites- 
imal condition in time (the jump differential equation 
'^t[^]p ~ [<^]p)' without simultaneously imposing the 
initial value of [vPJp that is consistent with the singular 
source [Eq. ([l6|]. Therefore, the prescription to elimi- 
nate Jost junk solutions is simple: To enforce the initial 
value [ ^ ]p (to) along the evolution so that the sources 
are correctly implemented. 

In light of this discussion the results and conclusions 
in U are correct. In particular, it is concluded there 
that Jost junk solutions are not numerical artifacts but 
rather they are related to the implementation of the sin- 
gular source term. However, the discussion about the 
specific reason underlying the appearance of a Jost solu- 
tion is not conclusive, as the failure to enforce the initial 



condition ( 16 ) is not identified as the underlying cause. 



Nevertheless, a way of avoiding the problem is proposed 
in [7] , consisting of a redefinition of ip |16) 



dt^ = dr^cj), (19) 

with J,j, = [vl'jp. No Jost solutions are found in the 
implementation of these equations with trivial ID for 
(^',(^, 0), and it is concluded in [7] that such kind of 
contamination is not generic. The reason is clear from 



the conclusions we have reached: system ( 19 ) implements 



the finite condition (|13|) on including (|16|), whereas 



system (17) does notT^ctually, system ( |19[ ) is equivalent 
to the second-order model Eq. ([T]), and this statement is 
independent of any discussion about ID [T7] . An alterna- 
tive way to understand the contrast between systems ( 17 ) 



and ( 19 1 is to note that trivial ID for ip and (p have (cf, 
[17]) completely different content. Whereas the ambi- 
guity stemming from the inconsistency between singular 
sources and trivial ID was resolved in the case of (17 1 by 



keeping trivial ID and ignoring the 6' term (by ignoring 
the value of ['^]p (to)), in system (19) one instead prior- 
itizes the correct implementation 01 the singular source 
from t = to at the (minimal) price of modifying the trivial 
ID on (p through ([l8| and 'p\t=ta = 0, which is equivalent 
to (16). In relation to this we make a statement that we 



(18) 



discuss later: Preserving the singular source at late times 
[and this, in certain implementations of ([!]) like (17) and 
(19), depends critically on the correctness of the source 
at t = to] systematically removes the Jost junk solution, 
independently of the ID. 

The problem discussed here is generic in the sense that 
it affects any linear system with distributional Dirac- 
delta sources in which (i) the sources are translated into 
jumps of the fields, and (ii) these jumps are implemented 
infinitesimally in time through an evolution equation. 
Consistency with the original system of equations de- 
mands the explicit enforcement of the initial value of the 
jump, otherwise a Jost junk solution will be present. 

To illustrate this point, we consider the computation of 
the self force in the simplified scenario corresponding to a 
charged scalar particle orbiting a nonrotating MBH. The 
spherical harmonic modes of the retarded scalar field, 
^im, satisfy the model Eq. ([!]) with the following iden- 
tifications: ^E* — r^im, F = 0, and G = f^^Sem [see, 
e.g. [TU], for the expressions of V{r) and Simit, r)]. With 
this, we can use the first-order reduction of Eqs. ^ or 
the one presented in [11] based on characteristic vari- 
ables, (^,M,w), with (u, w) = {(j> — (p, (j> + ip). Applying 
the splitting of the PwP approach [Eqs. ([9|-(11)] we ob- 
tain homogeneous equations and jump conditions: 

dt^^^{u^+v^)/2, [*]p(i) = 0, 
dtu^ ^~dr*u^ ~V{r)-^^ , [u]p{t) = -G{t,rp), 
dtv^ = dr-^v^ - F(r)*± , [v]p{t) = G{t,rp) . (20) 

The crucial point here is the implementation of jumps 
[u]p (t) and [v]p (t). In [TT] two different strategies have 
been adopted: (I) Finite jumps [u]p{t) and [w]p(t) are 
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FIG. 1: Jost solution built by subtracting solutions to system 
(20) in approaches (I) and (II) for jump implementation (cf. 
text), and using trivial ID in both cases. 



directly enforced in the evolution and trivial ID is used; 
(II) The time derivatives of the jumps, d[u]p/dt and 
d [ w ]p /dt, are imposed as extra evolution equations using 
the method of lines. Here, trivial ID cannot be employed 
and one must "impose initially the values of the jumps 
[...], since during the evolution the only input about them 
is the information on their derivatives" (cf. discussion 
in Sec. V of [Hj). Instead, appropriately modified ID is 
used [note the similarity to the discussion of system ( 19 )] . 

Therefore, a natural question is whether or not a Jost 
junk solution appears when taking the difference between 
the solutions obtained by implementing approaches (I) 
and (II) but using trivial ID in both cases. Figure [l] shows 
the results obtained from the numerical implementation 
with the techniques discussed in Ref. [TT|, which con- 
firm our conclusion that Jost junk solutions correspond to 
incorrect implementations of the distributional sources, 
rather than to (trivial) ID inconsistent with the sources. 



Hitherto, we have focused mainly on consistency is- 
sues at initial times. In order to get a deeper insight, 
we paraphrase the main discussion above in terms of so- 
lutions at late times. The only relevant element in the 
analysis at late times is the correct implementation of 
the sources, the choice of ID playing no role. In partic- 
ular, any intermediate-time data can be taken as valid 
ID. From this perspective, the crucial impact of the fail- 
ure in passing the initial value of the field jumps [TS] is 
the incorrect implementation of the sources, which differ 
from the correct ones at late times. Then, one is solving 
a different problem. 

The reason why a Jost solution appears when evalu- 



ating ^';J^st [Eq. ([4|] in [7^ is because the smooth time 
switch-on function, a(t,T), guarantees the correct imple- 
mentation of the late time sources in the solution 'I'smooth 
when adopting a strategy in which the source is im- 
plemented through an evolution equation for the jumps 
with zero initial values. To illustrate this, let us con- 



sider a field x with jump condition [x]^ = Jxi^) and its 
smoothed version [x]^ = a{t,T)J^{t). If we implement 
the jumps using evolution equations, d[x]p / dt = J'^ and 
d[x]p/dt — [oiJ-^y respectively, both with zero initial 
values, [x]p (to) = [x]p (*o) = 0, we obtain 

[X]p - Jx(i)- Jx(*o), 

[x\p = a{t,T)J^{t) - a{to,T)J^{to) 

= a{t, T)J^{t) ~ J^(t) (for i > r) . (21) 

This shows conclusively that 5'smooth solves Eq. ([T]) 
at late times with the correct source term, whereas 
^Impulsive is contaminated by a Jost solution [cf. point 
(ii) in the discussion after Eq. Q]. 

Time switch on in self-force calculations. Although 
the smooth switch on is not the critical element in the 
discussion of the Jost solution [19] . it certainly improves 
the consistency between the source and the ID [which 
can be critical in certain implementations of the model 
Eq. ([T])]. First, it reduces the initial burst of junk radia- 
tion, addressed in detail in [7]. Second, it improves the 
calculation of the self force in certain implementations. 
In [7; a similar point was addressed: The calculation of 
the gauge invariant gravitational-wave fluxes of energy 
and angular momentum. Regarding the self force itself, 
we have computed it for the case of a charged scalar par- 
ticle in circular orbits using the mode sum regularization 
scheme El. The retarded field is evolved using the sys- 



tem (20) and trivial ID in the approach (I), i.e. with the 
direct enforcement of the /imte jumps [u]^ and [v]^ (this 
is an example of a correct implementation of the source 
term while using trivial ID and, as expected, no Jost so- 
lution appears). On the other hand, the instantaneous 
switch on of the source term produces high-frequency nu- 
merical noise with a very slow time decay (the decay time 
scale is much larger than the orbital time scale) and dete- 
riorates significantly the accuracy of the self force. This 
noise can be eliminated by increasing the resolution of 
the computations or by using a numerical filter. How- 
ever, a more efficient and better adapted method is the 
smooth switch on of the source, which produces dramatic 
improvements in the computation of the self force. In 
Fig. [2] we illustrate the calculation of the self force, with 
and without time switch-on of the sources, for a parti- 
cle at the last stable circular orbit (r^ = 6i\f ) after two 
orbits, i.e. for t > 2 • To^b = 2 • 27r(r3/M)3 ^ 184.69M 
(the self force calculation is physically meaningful at late 
times, once the initial unphysical burst has been radi- 
ated away from the particle). Outgoing boundary condi- 
tions discussed in [TUl E] are imposed at the boundaries 
[?'*innor ~ -292M, r*outGr « 307Af], whose location guar- 
antees their causal disconnection from the particle up to 
the final evolution time tfin = 250M. 

To conclude, we have analyzed the risk of contamina- 
tion of time-domain self-force calculations by Jost junk 
solutions as a consequence of using trivial ID inconsistent 
with the distributional sources, which has produced some 
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FIG. 2: Radial and time components of the self force as 
a function of time, with ("On") and without ("Off") time 
switch-on of the sources (rp = 6M). Note that the Off case 
does not even allow the recovery of the second significant digit 
in $^ = 1.677 ■ 10"*, $t = 3.609 • lO-"* (see e.g. [TO]). 

confusion in the community. We have shown that no Jost 
solution is generated as long as the distributional Dirac- 
delta sources are faithfully implemented at late times. 
Indeed, the inconsistency between trivial ID and distri- 
butional sources introduces a genuine ambiguity in the 
evolution system: either one modifies the sources to make 
them compatible with the trivial ID, or one accepts the 
use of inconsistent ID and prioritizes the implementation 



of the correct sources (at late times). In this specific 
sense, the conclusion in [7] stating that the presence of 
Jost solutions is not a numerical artifact, but actually 
depends on the implementation of the system, is correct. 
However, it is not a generic feature of these systems and 
hence can be easily avoided. Here, we have first iden- 
tified unambiguously the specific origin of Jost solutions 
in implementations of jump conditions due to Dirac-delta 
sources, namely, the failure to enforce the correct initial 
value in the evolution equations for the jumps. Second, 
our final conclusion is that no contamination of the re- 
tarded solution by Jost junk solutions happens as long as 
(late time) sources are correctly implemented, and this is 
conceptually independent of the ID and/or the use of a 
smooth time switch-on. As a by-product of the analysis, 
we have shown that the use of a smooth time switch-on of 
the sources [71 [12] prevents the high-frequency noise that 
spoils self-force calculations in some numerical schemes. 
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